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A quantitative description is given for the densification process of glasses resulting from glass formation at 
elevated pressures. Phenomenological relations are derived, or justified, which allow estimation of the densifica- 
tion rate k' (with respect to formation pressure) from various thermodynamic quantities and glass transition 
behavior. In addition, the estimation of K l may be facilitated by the application of the hole theory of Simha and 
Somcynsky. Using these relations k' is estimated, and the results from the different methods are compared for data 
from 23 different organic polymers with glass transition temperatures ranging from 150 to 455 K. The amount of 
densification appears to be limited by the apparent convergence of the glass temperature and effective decomposi- 
tion temperature with increasing pressure. Some estimates of limiting values are presented. Finally, changes of 
refractive index resulting from densification are estimated from the observed, or predicted, densification rates. 
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1 . Introduction 

The density of a glass, as well as certain other properties, 
depend upon the thermodynamic history by which the glass is 
formed. For example, as shown schematically in figure la, 
an amorphous polymer subjected to an elevated pressure in \ 

the melt, followed by isobaric cooling at constant rate to a 
temperature well below the glass temperature, T g , and then 
depressurized, will have a larger density than that obtained 
by isobaric cooling at the same rate at atmospheric pressure 
to the same temperature in the glass. From the former proce- 
dure the pressure induced densification rate is defined as 



liquid; polymer; pressure; PVT; pyrolysis; 



k' = -a/v)(dv/dP') TJ >, k 



(i) 



where V is the volume at temperature T and pressure P, and 
P is the formation pressure maintained during constant rate 
of cooling k. Note that this definition parallels the usual one 
for the isothermal compressibility, 



K= -{l/V)(dV/dP) T j» Je , 



(2) 



the difference being that P andP' are interchanged. 

It is expected that the final depressurized volume in the 
glass will lie between the atmospheric and the pressurized 
value, as shown in figure la. It is then clear that the 
inequality, 
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Although we may intuitively expect this relation to hold, as 
apparent from experiment, we do not know of any proof. 

A well known alternative method of densifying glasses is 
simply to decrease the cooling rate as illustrated in figure lb. 
In principle one can obtain the same volume in the glass by 
this procedure as by elevating the pressure, except that the 
times required for the former are much longer. For example, 
it is estimated [l] 1 that a poly (vinyl acetate) glass obtained by 
isobaric cooling at 800 bar 2 in 8 hours would require 500 
years to reach the same volume at the same terminal tempera- 
ture by cooling slowly at atmospheric pressure. It should be 
recognized, however, that the states of glasses at, the same 
volume, temperature, and pressure, but obtained through 
different histories, are not necessarily the same. As pointed 
out by Bree and coworkers [2], volume changes during iso- 
baric-isothermal volume relaxation [3] have a large effect on 
relaxation times for creep compliance, whereas almost no 
effect is observed from volume changes obtained by isobaric 
cooling at elevated pressures. Accordingly, it appears that 
the state of a glass is not determined by its volume, tempera- 
ture, and pressure alone. Moreover, pressure induced densi- 
fication does have an influence on physical properties. Ac- 
cording to the data of Dale and Rogers [4] over a 5 kbar 
range, the compressive modulus of polystyrene appears to 
increase slightly with formation (or molding) pressure, level- 
ing off at higher pressures, with the yield stress going through 
a maximum between 1 and 1.5 kbar. Wetton and Money- 
penny [5] have studied the dynamic mechanical and dielec- 
tric properties of several polymeric glasses formed at pres- 
sures up to slightly beyond 5 kbar. Both the real part of 
Young's modulus and its loss tangent, as well as the real part 
of the dielectric constant, increase with formation pressure. 
For polyvinyl acetate) McKinney and Goldstein [1] have 
observed a 3 percent increase in the bulk modulus at ° C, 
corresponding to a formation pressure of 800 bar. This differ- 
ence increases with decreasing temperature. 

Thermal properties also seem to depend on the amount of 
pressure induced densification. Although the heat capacity 
C p is found to be independent [6] of formation pressure, the 
enthalpy H seems to vary significantly at formation pressures 
above a certain value. According to Price [7], very little 
change in the enthalpy of poly(methyl methacrylate) is ob- 
served up to about 800 bar, followed by a nearly constant rate 
of increase of about 0.015 cm 3 /g up to 3 kbar, their maximum 
value. For polystyrene [8], the data have been evaluated as 
A// = //(densified) — //(normal) first decreasing slightly and 
then going back to zero at about 800 bar, followed by an 
increase with nearly constant slope up to the maximum 
pressure. Weitz and Wunderlich [9] have also observed this 
behavior and interpreted it in terms of two opposing mecha- 
nisms arising from holes and rotational isomers. It is not 
clear, however, that the apparent negative values of AH 
obtained by experiment are significant. 

The purpose of this paper is to describe the thermody- 
namics of the pressure induced densification process by 
applying both phenomenological and molecular theory. Sim- 
ple phenomenological relationships are derived between the 
densification rate k' and other thermodynamic properties for 
which values are more readily available in the literature. 
Moreover, it is shown how the hole theory of Simha and 
Somcynsky [10] may be used to facilitate the estimation of the 



densification rates for polymers using a minimum amount of 
experimental information. In both cases the derived relation- 
ships are tested using appropriate experimental data. An 
example of the utility of these results is demonstrated by 
estimating the change in the index of refraction correspond- 
ing to changes in molding pressure, assuming that the index 
of refraction is related to the volume by the Lorentz-Lorenz 
equation. The results have potential application to the adjust- 
ment of the refractive indices of lenses by varying the mold- 
ing pressure. 

2. Phenomenological Relationships 

Two types of thermodynamic histories, shown schemati- 
cally in figure 2, are pertinent to the development of the 
phenomenological relationships for pressure induced densifi- 
cation. In the first (Fig. 2a) the PVT surface of the glass is 
obtained from repeated isobaric cooling runs at the same 
constant rate, but at difference pressures, with all pressure 
changes occurring in the melt prior to each run. This proce- 
dure is called the variable formation history because the 
structure of the glass is different for each experimental 
pressure (which is the formation pressure, since P — P'). The 
glass transition at each pressure is assumed to occur at a 
constant mean relaxation time. Hence, the intersection of the 
liquid and glass PVT surfaces gives the proper T g (P), from 
which dTg/dP is expected to approximate that obtained from 
the dynamic mechanical and dielectric frequency-tempera- 
ture-pressure superposition. On the other hand, as a conse- 
quence of the varied structure, the glass PVT surface is not 
proper in the thermodynamic sense. 
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FIGURE 2. Schematic illustration of two thermodynamic histories used to 
form glasses. 

(a) V ariable formation, (b) Constant formation. 



284 



With the other history (fig. 2b) the glass is formed also by 
isobaric cooling at constant rate at an arbitrary pressure, P' 
(which is usually atmospheric, but elevated in figure 2b to 
illustrate the more general case). At temperatures well below 
T g , where viscoelastic relaxation times are large in compari- 
son to effective experimental times, a thermodynamically 
reversible PVT surface for the glass is obtained by observing 
the volumetric response to "fast" changes in temperature and 
pressure. Since all of the data in the glass pertain to the same 
P', the PVT surface gives the proper values of the derivable 
thermodynamic quantities (for example, thermal expansion, 
isothermal compressibility, and internal pressure). The inter- 
section of the liquid and glass surfaces defines the fictive 
temperature Fj(P, P'). 

The principal distinction in procedure between the two 
histories is that with variable formation all pressure changes 
are made in the melt, whereas with constant formation they 
are made in the glass. Note that the number of independent 
variables is different for T g (P) and Tj(P, P'). The redun- 
dancy of using two arguments in the former arises from the 
fact that the formation and experimental pressure are always 
identical. Accordingly, T g may be regarded as a special case 
of 7 i when P = P' . The implicit argument k is deleted here 
because only one value appplies to these discussions for each 
case. For further details and interpretations of these histo- 
ries, see Ref. [ 1 ]. 

In all of the schematic diagrams in this paper the glass 
transition is shown as a discrete intersection. With isobaric 
cooling at constant rate through T g a gradual transition proc- 
ess is observed. The discrete intersections shown correspond 
to those obtained by extrapolation of the equilibrium isobars 
and the isochronal (nonrelaxing) ones for the glass. 

We now proceed to evaluate the thermodynamic diagram in 
figure 3, in order to determine relationships for K f in terms of 
other measured quantities. Volume A is obtained by isobaric 
cooling at constant rate and atmospheric pressure (P = P' = 
0). Volume/? is reached by pressurizing to P = P f = AP in 
the melt, followed by isobaric cooling at the same constant 
rate as for A, with subsequent depressurizing in the glass at 
the same terminal temperature as for A. Note that T g (in lieu 



P'=AP 




FIGURE 3. Schematic illustration of the procedure used to derive the 
densification equation [See eqs (3) and (4)], and the distinction between T g (P) 

and tMP,F). 



of Fj) applies here, since the transition is observed at the 
formation pressure (P = P f ) in both cases. The isobaric 
extension of V b with increasing T (see dashed line) to its 
intersection with the liquid line yields the fictive temperature 
T](P,P')=T](0,AP). 

In summing the thermodynamic contributions for small 
changes in T and P near T g (P) in the range where linear 
approximations are valid, we find 

V A = V (l - a l *T 1 - a g AT 2 ) 
V B = V {1 - k l AP - ai \AT t - (dT g /dP)AP] (3) 
- a g [AT 2 + (dT g /dP)AP] + Kg AP} 

where a is the usual isobaric thermal expansivity, the A's 
indicate differences as shown on figure 3, and the subscripts / 
and g pertain to liquid and glass. For small changes eq ( 1 ) 
may be written in the form 

«' = (V A - V B )/(VoAP) (4) 

Substitution of eqs (3) for V A and V B yields 

k' = Ak - AadTjdP (5a) 

where the A's here indicate the usual differences in the 
respective quantities between liquid and glass. Since 



drt/dP= Ak/Aol 



(6a) 



along Fj, [11], k' may be expressed in terms of the differ- 
ence between the two transition rates, i.e., 



k' = Aa(dT^/dP - dTjdP) 



(5b) 



Equation (5a) may also be written as an Ehrenfest-type 
relation, viz., 



dTg/dP = Ak/Acx - k'/A<x 



(5c) 



which is 
equality 



consistent with the experimentally observed in- 



dT a /dP < Ak/Aol 



provided the densification rate is non-negative. Expressions 
fully equivalent to eqs (5) have been derived by Goldstein 
[11] and given previously in Ref. [12]. From eqs (5) it is 
evident (as also pointed out by Goldstein [11]) that the 
necessary and sufficient condition (assuming Aa =/= 0) for the 
PVT surface to be independent of formation pressure is the 
validity of the first Ehrenfest equation 



dT g /dP = A/c/Aa. 



(6b) 



The analogous argument applies to the entropy surfaces. 
Since the second Ehrenfest equation, 



dTjdP = T a VAa/AC P 



(7) 



where Cp is the usual heat capacity at constant pressure, 
appears to be a good approximation [11, 13], there should be 
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a single entropy surface with respect to formation pressure in 
contrast to the manifold/surface observed for volume. This 
view is confirmed by the DSC 3 measurements of Yourtee and 
Cooper [6] on normal and densified polystyrene, which reveal 
no significant effect on the thermal properties of glasses by 
vitrification at elevated pressures. The authors did find some 
differences in the thermal behavior between these properties 
and those from vitrification by isothermal compression; how- 
ever, these were attributed to inhomogeneous freezing proc- 
esses during compression. Accordingly, if eq (7) is a good 
approximation, it leads to a convenient experimental deter- 
mination of the initial (P — 0) value of dT JdP through 
volume-temperature and heat capacity measurements re- 
quired at atmospheric pressure only. Equation (7) will be 
tested by means of experimental data later in this paper. 

As stated above k' may be determined (near 7^) from the 
values of Aa, dTJdP, and Ak using eq (5a). The relative 
difficulty in obtaining these quantities experimentally in- 
creases in the order given above, as does the difficulty of 
obtaining their values from the literature. For these reasons it 
is desirable to be able to estimate Ak (or A/c/Aa) independ- 
ently of existing PVT data. It will be shown how the hole 
theory of Simha-Somcynsky [10] may be used to arrive at 
values of dT] / dP = Ak/A«. 

As indicated previously, eqs (5) are based on several 
linearizations. It is assumed that the coefficients (Xi, a g , Kj 
and Kg are independent of pressure and temperature and that 
T g is a linear function of pressure. Thus strictly, the reference 
temperature in the glass as well as the initial temperature in 
the melt should be appropriately close to T g . Moreover, the 
pressure P' should be appropriately small. In the Appendix 
the general relationships are developed, based on the equa- 
tions of state of the liquid and both glasses. 

As an example, integral relations are evaluated over the 
two paths shown on figure la for PVAc, for which extensive 
data are available [1], and the Tait parameters are known [14] 
for the liquid and both glasses. The results are tabulated and 
compared with the corresponding linear approximations. 

3. Application of Molecular Theory 

Th$ hole theory, which is used here to estimate the values 
of dT]/dP, is a corresponding states theory based on a lattice 
model. The partition function is defined in terms of a single 
ordering parameter, the hole fraction h, which gives the ratio 
of the number of vacant to total sites, each of which may be 
occupied by a polymer segment. The corresponding states are 
given in terms of the reduced (universal) variables 



and the equilibrium constraint (dZ /dh) T} v = 0, the following 
equilibrium equations [10] are obtained, respectively: 



T = TIT*, P = PjP*, V = V/V* 



(8) 



where 7 1 *, P*, and V* are the scaling factors applicable to 
each polymer. Although these are defined explicitly by the 
theory, they are usually derived from a superposition of 
equilibrium PVT data along the master curves evaluated from 
the theory. For an illustration of this procedure, see Ref. 
[14]. 

The partition function Z is expressed uniquely in terms of 
the three independent variables T, V, and h. From the 
thermodynamic definition 

P = -kT[d€nZ(T, V, h)/dV] T 



(9) 



py/f = [1 - 2- 1/6 r (r^)" 1/3 ] _1 

+ (2y/t)(yV)- 2 [l.0ll(jV)- 2 - 1.2045] 
(s/3c)[(s - 1)A + y-Hn{\ - y)] 

= [2- 1/6 r (r^)" 1/3 - i/3][i - 2- 1/6 r W- 1/3 ]- 1 (10) 

+ [y/(6T)](yV)- 2 [2A09 - 3.003( r K)" 2 ] 

where y = 1 — h is the fraction of occupied sites, and s and 
3c are the number of segments per molecule and the external 
degrees of freedom per molecule, respectively. As in pre- 
vious work, we take s/3c = 1. Note that the term (s — l)/s in 
eq (10) approaches unity for large molecules. 

A basic assumption sufficient for the application of the 
hole theory to our densification model is 



dTJ/dP = (dT/dP) h . 



(11) 



3 Differential scanning calorimetry 



Gee [15] has shown that such an equation is valid for a single 
ordering parameter which is frozen in the glass. However, 
since h has been found to vary slightly with temperature and 
pressure in the glass [16, 17, 18], eq (11) must be revaluated 
to assess its validity for the more general case. 

Consider the single-valued function V = V (T, P, h) for 
which, by the usual definitions, 

-k = {d€nV/dP) T 

= (d€nV/dP) Tth + {d£nV/dh) T , P {dh/dP) T 

a = (d£nV/dT) P 

= (d€nV/dT) P>h + {d£nV/dh) T , P {dh/dT) P . 

Since there are three independent variables (in the general 
case, the derivatives with two fixed arguments (suscripts) 
are the same for liquid and glass. (For the glass it is under- 
stood that these derivatives pertain to constant P' and A:.) The 
differences become 

-Ak = (dfnV/dh) T , P [(dh/dP) T , € - (dh/dP) T J 

Aa = (d€nV/dh) TJ >[(dh/dT) P j - (dh/dT) P , g l 

where the subscripts £ and g again pertain to liquid and 
glass. Recalling that d7> /dP = Ak/Aq:, the ratio of the 
above equations is 

dTpdP = (dT/dP) h , ( F P /F T (12) 

where the "freezing fractions" F p and Ft are 
F P = 1 - (dh/dP) T J(dh/dP) T , ( 
F T = 1 - (dh/dT) Pu ,l(dh/dT) P , ( (13) 

as defined in ref. [14]. Note that whenF^ = F P , eqs (11) and 
(12) coincide. 
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Since the variables T, P, 
follows that 



and h are continuous at T t 

g ' 



(dT/dP) h j = (dT/dP) hiB = (dT/dP) h 
along this transition line. From partial differential equations, 



(dT/dP) h 



■(dh/dP) T j/(dh/dT) P , e 
- (dh/dP) T J(d/h/3T) Tjl , 



it follows that 



(dh/dP) T J(dh/dP) T j = {dh/dT) P J{dh/dT) P j, 

which, as seen from eq (13), is tantamount io F T = F p . 
Accordingly, since eqs (12) and (13) coincide, the validity of 
eq (12) is extended to a single ordering parameter which need 
not be "frozen" in the glass. 

To our knowledge both of the above freezing fractions have 
been evaluated for only two systems, namely polyvinyl ace- 
tate) [16] and selenium [19]. According to the best analysis 
given in ref. [16], F P = 0.88 and F T = 0.82 for which the 
ratio F p/F V = 1.07, which corresponds to a 7 percent 
discrepancy in eq (11) for polyvinyl acetate). A similar 
conclusion follows for Se. Since the above analysis shows that 
F t — F p, these differences are taken to be artifacts resulting 
from numerical inaccuracies. 

The next step is the evaluation of (dT/dP)^ at equilibrium. 
From simultaneous numerical solutions of eqs (9) and (10), 
values of h = 1 — y are obtained at a given set of reduced 
temperatures and pressures. For computational purposes it is 
convenient to replace (3T/dP) h by the ratio —(dh/dP)f/ 
(dh/3T)p. With constant increments Ax(x = T or P), it is 
easily shown for a quadratic dependence of y on x that 

(dy/dx)i = (y i+l - y i . 1 )/(2Ax). 

This procedure is used to generate a set of (dh/dP)f and 
(dh/dT)p values over the desired range from the sets of 
quadratic arcs defined by three adjacent points. Using a least 
squares fit, the approximation 

(dt/dP) h = 0.00502 + 0.198 T+ 31 .476r 2 (14) 

is found to be accurate within a residual standard deviation of 
0.2 percent at atmospheric pressure over the range 0.01 < T 
< 0.04. From eqs (8) and (11) 



dfl/dP = T*(dt/dP) h /P*. 



(15) 



Substitution of eq (15) into eq (5b) gives the desired relation 
for the densification rate, 



k' = ba[T*(dt/dP) h /P* - dT g /dPl 



(16) 



where all quantities are evaluated at T = T g . 

Equation (16) may be rewritten (see for example eq (14) in 
Ref. [14]) as 



k' = -Aa[(dT/dh) P X dh/dP] 



(16') 



where the total derivative on the right hand side is to be taken 
along the T g (P) line. Provided the pressure coefficient of T g 
has been determined with sufficient accuracy, there appar- 
ently is no numerical advantage in using eq (16'). 

It is mentioned in the last section that eq (7) appears to be 
a good approximation for most polymers. Assuming this 
relation, we may estimate k' from volume-temperature and 
heat capacity data using the relation 

*c' = Aa[T*(df/dP) h /P* - T g V Aa/A6>|. (17) 

Olabisi and Simha [17] have shown for most polymers studied 
by them that the scaling factor P* may be determined from 
the other two by means of the empirical relation 

p* = (T*/V*) exp (1.319 - 1.493 X 10~ 4 r*) (18) 

where the dimensions are K, bar, and cm 3 /g. Thus it appears 
feasible to estimate k' from appropriate data at atmospheric 
pressure only. This possibility is tested later in this paper. 

According to Wunderlich [20] it is possible to estimate 
ACpat T g to within about ±2J/(mol — K) by applying the "rule 
of constant ACp". The molecular repeat units are broken up 
into fundamental units or "beads" which loosen up in the T g 
process. Each bead is assigned the value 11.3 J/(mol— K). 
The contributions of the beads to AC p are assumed to be 
additive. Accepting the validity of this rule, it appears possi- 
ble to obtain a crude estimate of k' from hole theory applied 
to volume-temperature measurements alone. 



4. Results 



4.1. Data Sources 

Although we refer usually to the original sources, there are 
collections of data on the pertinent quantities in the litera- 
ture, which are sometimes cited here. Extensive lists of 
polymers and their values of T g appear in refs. [21-23], 
where the last is restricted to fluorine-containing systems. 
Tables of T g and Act are included in Refs. [24-26], and T g 
and AC P in Refs. [27] and [28]. Reference [29], which is 
occasionally cited here, contains a more critical evaluation of 
C p data on polymers for which the values on the same 
substance are often based on averages from different sources 
over wide ranges of temperature. Lists of polymers and their 
scaling factors based on the hole theory appear in Refs. [17] 
and [30]; however, P* is not available in the latter. Pyrolysis 
data on polymers are contained in refs. [31-34]. These are 
useful to prevent degradation during the densification process 
and to optimize the amount of densification. Finally, ref. [35] 
gives an extensive list of refractive indices for polymers. 

The number of digits for the values given in the subsequent 
tables is not intended to be an indication of precision or 
accuracy. Usually these numbers correspond to those given by 
the data sources. It is our opinion that most of the entries in 
these tables have more digits than can be justified as signifi- 
cant. 

Table 1 gives the lists of polymers studied, abbreviations 
used here, and their values of T g . In all tables the sequence is 
in order of increasing T g . 
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Table 1. List of polymers studied, abbreviations, and glass temperatures 



Abbreviations 


Polymer of 


TVK) 


PDMSi 


Dimethyl siloxane 


150 


PIB 


Isobutylene 


202 


NR 


Natural Rubber 


204 


PP 


Propylene 


244 


SBR a 


Styrene-butadiene 


257 


PMA 


Methyl acrylate 


282 


PnBMA 


rc-butyl methacrylate 


293 


PVAc 


Vinyl acetate 


304 


i-PMMA 


Isotatic methyl methacrylate 


320 


PEMA 


Ethyl metacrylate 


337 


PET 


Ethylene terephthalate 


340 


PVC 


Vinyl chloride 


349 


P4POS 


4-phenoxystyrene 


352 


P3CS 


3-chlorostyrene 


363 


PS 


Styrene 


363 


a-PMMA 


Atatic methyl methacrylate 


378 


PCHMA 


Cyclohexyl methacrylate 


380 


P4M0S 


4-methoxy styrene 


381 


P4CS 


4-chlorostyrene 


400 


PoMS 


o -methyl styrene 


404 


Pcarb 


Carbonate of bis(phenol A) 


416 


PaMS (1) 


Q:-methyl styrene (67% syndiotatic) 


448 


PaMS (2) 


a-methyl styrene (95% syndiotatic) 


455 



55 percent Styrene. 



4.2. Scaling Factors 

Table 2 gives the scaling factors based on the hole theory 
of Simha and Somcynsky [10]. These are determined through 
superposition of experimental equilibrium data on each poly- 
mer with respect to the theoretical equation of state. In this 







Table 2. Polymer scaling factors 










j* 


v* 


p* 


P*calc 


AP*/P 


Polymer 


Ret. 


K 


cm 3 /g 


bar 


bar 


% 


PDMSi 


36, 


30 a 


7893 


0.9602 


5061 


9461 


87 


PDMSi 




37 


7893 b 


.9602 b 


5061 


9461 


87 


PIB 


38, 


30 a 


11220 


1.0902 


— 


7208 


— 


PIB 




37 


11220 b 


1.0902 b 


7316 


7208 


-1.4 


NR 




39 


8344 


1.0353 


— 


8672 


— 


PP 




40 


8966 


1.1230 


8437 


7829 


-7.2 


SBR C 




41 


9800 


0.9892 


— 


8577 


— 


PMA 


42, 


17 


9200 


.7925 


— 


10990 


— 


PnBMA 


43, 


17 


9988 


.9299 


8456 


9402 


11 


PVAc 


1, 


14 


9419 


.8141 


9380 


10600 


13 


PVAc 




37 


9419 b 


.8141 b 


9158 


10600 


16 


i-PMMA 




44 


11170 


.8160 


10090 


9659 


-4.3 


PEMA 


45, 


30 a 


11540 


.8862 


— 


8694 


— 


PET 




46 


10870 


.7406 


— 


10830 


— 


PVC 


47 d , 


17 


11320 


.7105 


10350 


10990 


6.2 


PVC 




48 


11363 


.7083 


9783 


11000 


12 


PS 




15 


12880 


.9601 


6628 


7333 


11 


PS 


47 d , 


17 


12700 


.9625 


7638 


7409 


-2.9 


PS 




49 


11630 


.9480 


7881 


8082 


2.6 


PS 




50 


12680 


.9598 


7453 


7440 


-0.2 


a-PMMA 




48 


11490 


.8220 


8987 


9403 


4.6 


a-PMMA 


43, 


17 


11920 


.8370 


9147 


8984 


-1.8 


a-PMMA 


47 d , 


17 


11890 


.8350 


9303 


9024 


-3.0 


PCHMA 


43, 


17 


11290 


0.8906 


8382 


8786 


4.8 


PoMS 




50 


12740 


.9762 


7458 


7285 


-2.3 


Pcarb 


51, 


30 a 


12130 


.8100 


— 


9156 


— 


PaMS (1) 


52, 


30 a 


12700 


.9152 


— 


7792 


— 


PaMS (2) 


52, 


30 a 


12790 


.8433 


- 


8403 


- 



a T* and V* only are determined in this reference. 

b Volume-temperature data at atmospheric pressure taken from previous 
listing. 

c 55 per cent Styrene. 

d For additional comments of interpretation on these experimental data 
and evaluation of Tait parameters, see Ref. [53]. 



work the scaling factors are used solely to estimate dT\ / dP 
= A/c/Aa for each polymer using eq (15). When two numbers 
appear in the reference column (in table 2 only), the first 
applies to the data source, and the second to the work by 
which the scaling factors are evaluated. When only one 
number appears, the scaling factors are either evaluated in 
the reference given, or by us. 

Two values of P* for each polymer (or row) usually appear. 
The first of these (P*) is determined in the usual way through 
superposition as mentioned above. The second (P*aic) is 
obtained from eq (18). When volume-temperature data are 
available at atmospheric pressure only, it is necessary to use 
eq (18) to estimate P*. With the exception of PDMSi, P* and 
^caic agree to within 17 percent with an 8 percent relative 
standard deviation of differences over 17 pairs. With PDMSi 
the disparity of 87 percent is outstanding, and it is to be 
noted that the reduced glass temperature lies significantly 
outside the range for which eq (18) was deduced. Similarly, 
V* and T* are obtained at considerably higher temperatures 
than those employed here. A decrease of T* by 7.2 percent 
and a concomitant decrease of K* by 2.4 percent over 100 K 
has been estimated [30] for this polymer. Accordingly, the 
scaling factors cannot be assigned significant constant values 
over the experimental range. 

In order to obtain some measure of the uncertainty in the 
scaling factors, several data sources on each polymer are 
sometimes included. 

4.3. Densification Rates from PIT Data 

Table 3 gives the results of calculations of the densification 
rates from PVT data without recourse to molecular theory. k\ 
is determined from the definition: 



«\ = (v A - V B )/V A P' 



(4') 



which is identical to eq (4) setting AP = P\ except that V A 
replaced V . The difference between the values of k[ deter- 
mined from eqs (4) and (4') are insignificant in comparison 
with experimental uncertainty. k' 2 is determined from eq 
(5b). Note that there are only two polymers, PVAc and PaMS 
for which we found sufficient information to determine both 
k[ and Ko. Although the two methods are not necessarily 
fully equivalent because of the assumptions used to derive 
eqs (5), the agreement in both cases is good. In instances of 
more than one set of values per polymer, it is clear that the 
deviations in dTj/dP have the largest effect on the uncer- 
tainty of K f 2 . These apparent discrepancies are usually con- 
sistent with the differences in A/c. With polystyrene the 
maximum deviation in dT\ /dP is 38 percent compared with 
those for dT 9 /dP, 23 percent and Aa, 10 percent. Since k' 2 
involves the difference between the two transition rates, its 
maximum deviation is magnified to 56 percent with its rela- 
tive standard deviation over the five values being 26 percent. 
It is interesting to note that the direct method giving k[ which 
one might expect to be more reliable, yields values for which 
the maximum deviation (for polystyrene) is 86 percent with 
relative standard deviation over seven values being 35 per- 
cent. The ratio of average values, Kjj/Jci, is 1.6. These 
discrepancies are a measure of the difficulties in obtaining 
reliable PVT data on glassy polymers. 

In many instances the values of Ac* at T Q and the required 
transition rates are not tabulated in the data sources and, 
therefore, had to be evaluated. The accuracy of these evalua- 
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tions may be considerably limited when the data are pre- 
sented in graphical form only. In ref. [47] the values of Ak 
are determined by a different definition from the one used by 
us. In our definition K g is taken to be an isochronal (non- 
relaxing) function of temperature and pressure and therefore 
must be derived from data at temperatures below (or pres- 
sures above) the glass transition region. Ak at T g is then 
obtained by extrapolation. This, apparently, was not done by 
Hellwege et al. [47], at least over the appropriate temperature 
range for the data to be effectively isochronal. The distinction 
between the two definitions is clearly illustrated by Boyer 
[57]. Our larger values of Ak are determined from the Tait 
parameters given in ref. [53], which apply to the data of 
Hellwege et al. [47]. Note that the values of dTt/dP from 
reevaluating [53] their data are in good agreement with most 
of the others on the same polymers except for PMMA. This 
discrepancy would be increased by using their value of Ak 
given in ref. [47], along with poorer overall agreement with 
the other two polymers. 

In ref. [1] the transition rates are given as tangent values 
along the transition lines at each experimental pressure. 
Here, we use the secant values dT g /dP and dTt/dP between 
and 800 bar. This procedure gives average values and is 
more consistent with other treatments. 

The smallest value of k' =0.7 Mbar -1 in table 3 applies 
to PnBMA. Such a small value implies that the first Ehrenfest 
equation is a good approximation for this polymer. [See eqs 
(5c) and (6b)]. 



Table 3. Densification rate calculations at atmospheric pressure from PVT 

data 



Polymer 


Ref. 


Ac* x 10 4 
K- 1 


dTJdP 
K/khar 


dT g VdP 

K/khar 


K> \ 
Mbar" 1 


k' 2 
Mbar" 1 


PnBMA 


43 


1.69 


20.4 


24.3 


_ 


0.7 


PVAc 


1 


4.32 


21.2 


42.4 


8.8 


9.2 


i-PMMA 


44 


3.49 


21.1 


35.2 


— 


4.9 


PVC 


47 


2.93 


13.5 


35. 5 a 


— 


6.4 


PVC 


48 


3.71 


14 


46 


— 


12 


PVC 


2 


— 


— 


— 


4.4 


— 


P4P0S 


5 


— 


— 


— 


7.2 


— 


P3CS 


5 


— 


— 


— 


5.8 


— 


PS 


15 


2.84 


31 


71 


— 


11 


PS 


47 


3.05 


30.0 


56. 4 a 


— 


8.1 


PS 


49 


3.15 


25.0 


49.8 


— 


7.8 


PS 


50 


2.97 


31.6 


74.2 


— 


13 


PS 


2 


— 


— 


— 


5.5 


— 


PS 


4 


— 


— 


— 


10 


— 


PS 


5 


— 


— 


— 


7.2 


— 


PS 


6 


— 


— 


— 


4.3 


— 


PS 


8 


— 


— 


— 


5.6 


— 


PS 


9 


— 


— 


— 


9.5 


— 


PS 


54 


— 


— 


— 


4.5 


— 


PS 


55 


3.12 


32 


73 


— 


13.8 


a-PMMA 


48 


3.1 


18 


71 


— 


15 


a-PMMA 


47 


2.95 


23 


35. 9 a 


— 


3.8 


a-PMMA 


43 


2.35 


23.6 


54.5 


— 


7.3 


a-PMMA 


2 


— 


— 


— 


5.9 


— 


a-PMMA 


54 


— 


— 


— 


4.9 


— 


PCHMA 


43 


3.38 


22.4 


59.8 


— 


13 


P4MOS 


5 


— 


— 


— 


7.1 


_ 


P4CS 


5 


— 


— 


— 


9.5 


— 


PoMS 


50 


2.71 


34.2 


73.0 


— 


11 


Pcarb 


2 


— 


— 


— 


5.1 


— 


PaMS b 


56 


3.98 


39 


52 


4.9 


5.2 



a The values of Ak here are not consistent with those given in Ref. [47], 
For distinction, see text and/or Ref. [57]. 

b Tacticity not known to correspond to listings in table 1. 



4.4. Theoretical Estimation of the Transition Rate 

dTf/dP 

As stated in section 3 the transition rate dTi/dP = Ak/ 
Aa applicable to the constant formation history, may be 
estimated from the value of (dT/dP) n at T g . Table 4 summa- 
rizes the results of these calculations. In all cases except for 
PDMSi and PIB the data encompass T g . With these two 
polymers the first reference for each applies to the source of 
data at atmospheric pressure, and the second, at elevated 
pressures. 

Equation (14) and the scaling temperatures and pressures 
as applied to eq (15) provide the requisite information. The 
distinction between the values of (dT/dP) n in columns A and 
B is that they correspond to P* and P* a i c , respectively, in 
table 2. Values of dTt/dP are included for comparison with 
those of (ST/dP)h in cases where there is sufficient experi- 
mental information. 

With natural rubber (NR) Ak is determined from dynamic 
compressibility data [60]. This involves the measurement of 
the adiabatic compressibility in a hydrostatic stress field 
alternating at low audio frequencies. The low- and high 
frequency limiting compressibilities are converted from adi- 
abatic to isothermal conditions, providing the difference Ak. 
The fact that Ak is determined at about 20 K above normal T g 
is expected to have no appreciable effect. 

Except for PDMSi, the corresponding values of {()T/()P) n 
are nearly the same in columns A and B. The discrepancy for 
PDMSi is explained by the factors mentioned earlier. (See 
sec. 4.2.). Excluding this polymer, the relative standard 
deviation of the differences between corresponding values in 
these columns is 6 percent, which is considered to be good 
agreement. For comparison with experiment, the residual 
standard deviation between corresponding values of {clT/dP) n 
(column A) and dn / dP is 18 percent. Poor agreement is 
noted for PnMBa, and two samples of a-PMMA. 

Over a single substance, for example polystyrene where we 
have four sets of values, the relative standard deviations with 
respect to the averages for {dT/dP) h (column A) and dTj / dP 
are 6 and 15 percent, respectively. That for the differences 
between corresponding values of these quantities is 11 per- 
cent. Thus, based on these simple statistics, the most serious 
limitation is not the inadequacy of the theory, but the uncer- 
tainty in the experimental determination of aT' /dP. The 
agreement between values of dT\ / dP from different investi- 
gators is even worse for a-PMMA. 

4.5. Application of Hole Theory 

After determining the values of (dT/dP) n for each sub- 
stance, k' may be estimated by eq (16), where dTJdP may be 
determined by means of PVT data, dynamic measurements at 
elevated pressures, or heat capacity and thermal expansion 
data, both at atmospheric pressure. 

PVT Data 

Since PVT measurements are often made by the variable 
formation history only (for example, polypropylene, ref. [40]), 
there is insufficient information to determine Ak, and hence 
dT\ jdP \ to be applied to eq (5b). Accordingly, this quantity 
is replaced by T*(dT/ dP) JP* leading to eq (16). The results 
of these estimates are given in table 5, where k' corresponds 
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to k' 2 in table 3. In all cases (dT/dP)^ is taken from column A 
of table 4. Since the expression for k' involves the difference 
between two transition rates, its value is very sensitive to (dT/ 
dP) h . This effect is reflected in the large standard deviation, 
2.5 Mbar -1 , with a relative value of 28 percent, for the 
differences over 14 pairs of corresponding values of k' 2 and 
k' in tables 3 and 5. The values of PaMS are not included in 
this calculation because there are insufficient data in ref. [56] 
to determine the scaling factors applicable to this particular 
sample. Based on the fact that the standard deviation of k' 2 
for polystyrene (table 3) over five values is 2.8 Mbar -1 
corresponding to 26 percent, the overall 28 percent value 
above appears to be dominated by experimental uncertainty. 



Table 5. Densification rate calculations using PVT data and hole theory 



Table 4 


Transition 


rates calculated from hole theory 






Ref. £ 


t 9 
x 10 5 


(of/dP) h 
x 10 5 


A 


B 


dT-tJ 


Polymer 


(oT/dP) h 
K/kbar 


(oT/dPh 
K/kbar 


dP 
K/kbar 


PDMSi 


36, 37 


150 a 


1900 


2014 


31.4 


16.8 


— 


PIB 


37, 38 


202 b 


1800 


1878 


28.8 


29.2 


— 


NR 


39 


204 


2445 


2868 


— 


27.6 


24 c 


PP 


40 


244 


2721 


3371 


35.8 


38.6 


— 


SBR 


41 


256 


2622 


3185 


— 


36.4 


— 


PMA 


42 


281 


3054 


4042 


— 


33.8 


— 


PnBMA 


43 


293 


2934 


3792 


44.8 


40.3 


24.3 


PVAc 


1 


304 


3226 


4416 


44.3 


39.2 


42.4 


t-PMMA 


44 


320 


2865 


3653 


40.4 


42.2 


35.2 


PEMA 


45 


338 


2928 


3780 


— 


50.2 


— 


PET 


46 


340 


3128 


4201 


— 


42.2 


— 


PVC 


47 


349 


3083 


4104 


44.9 


42.3 


41 


PVC 


48 


349 


3071 


4079 


47.4 


42.1 


46 


PS 


15 


362 


2811 


3546 


68.9 


62.3 


71.0 


PS 


47 


362 


2850 


3623 


60.2 


62.1 


56 


PS 


49 


365 


3138 


4223 


62.3 


60.8 


55.5 


PS 


50 


365 


2879 


3681 


62.6 


62.7 


74.2 


a-PMMA 


48 


376 


3272 


4520 


57.8 


55.2 


71 


fl-PMMA 


47 


378 


3179 


4312 


55.1 


56.8 


35.9 


a-PMMA 


43 


378 


3171 


4295 


56.0 


57.0 


54.5 


PCHMA 


43 


380 


3366 


4735 


63.8 


60.8 


59.8 


PoMS 


50 


404 


3171 


4295 


73.4 


75.1 


73.0 


Pcarb 


51 


423 


3487 


5020 


_ 


66.5 


— 


PaMS d 


56 


440 


3465 


4967 


— 


81.0 


52 


PaMS (1) 


52 


448 


3528 


5118 


— 


83.4 


— 


PaMS (2) 


52 


455 


3557 


5189 


- 


86.5 


- 



a T g taken from ref. [58]. 
b T g taken from ref. [59]. 

c A/c determined from dynamic compressibility data [60] on vulcanized 
natural rubber with 12 percent combined sulfur. See text. 
d Scaling factors taken the same as for PaMS (1). 



Dynamic Data 

From the assumption that the value of dTJdP, approxi- 
mates that of (dT/dP)^, where (o is the angular frequency, 
frequency-temperature-pressure superposition of dynamic 
data, including dielectric and ultrasonic, may be used to 
determine dTJdP appearing in eq (16). The validity of this 
assumption is, of course, subject to the condition that (dT/ 
dP)^ for the T g -process is essentially independent of fre- 
quency. (Numerical comparisons between different experi- 
mental transition rates are made below.) 

The results of these calculations are given in table 6. Since 
values of (dT/dP) n in column A of table 4 involve fewer 
assumptions than those in column B, the former are used 
where there is a choice. The standard deviation of the differ- 
ences of k' over eight pairs, where there are values on the 



Polymer 


Ref. 


Aa x 10 4 
K- 1 


dTJdP 
K/kbar 


(oT/dP) h 
K/kbar 


k' Mbar ! 


PP 


40 


4.41 


20 


35.8 


7.0 


PnBMA 


43 


1.69 


20.4 


44.8 


4.1 


PVAc 


1 


4.32 


21.2 


44.3 


10 


i-PMMA 


44 


3.49 


21.1 


40.4 


6.7 


PVC 


47 


2.93 


13.5 


44.9 


9.2 


PVC 


48 


3.71 


14 


47.4 


12 


PS 


15 


2.84 


31 


68.9 


11 


PS 


47 


3.05 


30.0 


60.2 


9.2 


PS 


49 


3.15 


25.0 


62.3 


12 


PS 


50 


2.97 


31.6 


62.6 


12 


a-PMMA 


48 


3.1 


18 


57.8 


12 


a-PMMA 


47 


2.95 


23 


55.1 


9.5 


a-PMMA 


43 


2.35 


23.6 


56 


7.6 


PCHMA 


43 


3.38 


22.4 


63.8 


14 


PoMS 


50 


2.71 


34.2 


73.4 


11 


PaMS 3 


56 


3.98 


39 


81 


17 



Scaling factors taken for PaMS (1), Ref. [52]. 



same substances in table 3 including k\ for Pcarb, is 2.2 
Mbar -1 or 30 percent, which is about the same as the 
experimental uncertainty given above (2.8 Mbar -1 or 26 
percent) for polystyrene. This value is also about the same as 
the 28 percent value given for the PVT data even though data 
on different substances are involved. It is possible, however, 
that dTg/dP values determined from dynamic data, in partic- 
ular dielectric, where high resolution is obtained, are more 
reliable than PVT values. T g determinations from PVT data 
usually involve extrapolations which are not used in the 
superposition of dynamic data. 

Heat Capacity and Thermal Expansion 

The estimation of k' from heat capacity and thermal ex- 
pansion data is based on the apparent validity of the second 
Ehrenfest equation [eq (7)]. (This relation is tested in the next 
section.) The results of the calculations based on eq (17) are 
summarized in table 7. In this case the statistics may not be 
meaningful because there are only five values of k' which 
correspond to those in table 3 including K^ for Pcarb. PaMS 
is excluded for the reason given above. The standard devia- 
tion of the differences is 1.4 Mbar -1 or 17 percent, which is 
somewhat less than the experimental uncertainty (26 percent) 
based on polysytrene data (table 3). In view of the high 
experimental uncertainty for all methods, this method of 
estimating k' appears to be reliable, except for PDMSi and 
PaMS. 

In ref. [12] a negative value of k' for PDMSi [based on eq 
(17)] is reported. This is a surprising, but not necessarily an 
incorrect result. The analysis of this polymer is hindered by 
the lack of good thermal expansion data through T g , largely a 
consequence of the low temperatures required, and the strong 
tendency for this polymer to crystallize. The negative value of 
k' is obtained by using the value of Act = 10.28 K -1 from 
table I of ref. [26]. This value is based on the linear thermal 
expansion data of Weir, Leser, and Wood [58]. After a 
thorough examination of their results and consultation with 
Dr. Wood, it was decided that the temperature range for 
which V OLi was evaluated is too small and too remote from T g 
to evaluate Aa at T g . In order to obtain what we consider to 
be the best available estimate of Aa, we used the value of 
Va g - 2.7 X 10 -4 cm 3 /(g-K) from ref. [58], and a t = 8.7 X 
10 -4 K -1 and V g = 0.904 cm 3 /g from the density-temperature 
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Table 6. Densification rate calculations from dynamic data at elevated pressures and hole theory 



Polymer 


Ref. 


Aa x 10 4 K- 1 


(8T/dP) h K/kbar 


Ref. 


Method 


(oT/dPL K/kbar 


k' Mbar ! 


PIB 


37, 38, 59 


4.2 


28.8 


61 


U 


25 


1.6 


NR 


39 


5.40 


27.6 


60 


C 


24 


1.9 


PMA 


42 


3.7 


33.8 


62 


1) 


18 


5.8 


PnBMA 


43 


1.69 


40.3 


63 


D 


16.7 


4.0 


PVAc 


1 


4.32 


39.2 


13 


D 


22 


7.4 


PVAc a 


1 


4.32 


39.2 


64 


C 


20 


8.3 


PEMA 


45 


2.95 


50.2 


65 


D 


20 


8.9 


PVC 


47 


2.93 


44.9 


66 


D 


18 


7.9 


PVC 


47 


2.93 


44.9 


67 


S 


16.5 


8.3 


PS 


50 


2.97 


62.6 


68 


D 


32 


9.1 


a-PMMA 


43 


2.35 


56.0 


67 


S 


24.5 


7.4 


Pcarb 


51 


2.81 


66.5 


13 


D 


44 


6.3 



a All quantities derived for same sample. 

U Ultrasonic 

D Dielectric 

C Dynamic Compressibility 

S Dynamic Shear 



Table 7. Densification rate calculations from thermal expansion and heat capacity data and hole theory 



Polymer 


Ref. 


T g K 


V cm : 7s 


Aa x 10 4 K" 1 


(oT/dP) h K/kbar 


Ref. 


AC p J/(g-K) 


dT g /dP a K/kbar 


k' Mbar ' 


PDMSi b 


58, 30 


150 


0.904 


6.0 


31.4 


27 


0.30 


27.0 


2.6 


PIB 


59 


202 


1.072 


4.2 


28.8 


69 


.377 


24.1 


2.0 


NR 


39 


204 


1.023 


5.40 


27.6 


70 


.455 


24.8 


1.5 


ppc 


40 


244 


1.127 


4.41 


35.8 


71 


.51 


23.8 


5.3 


SBR d 


41 


256 


0.987 


3.92 


36.4 


72 


.456 


21.7 


5.8 


PVAc e 


1 


304 


.843 


4.32 


44.3 


16 e 


.50 


22.1 


9.6 


PVC 


47 


349 


.729 


2.93 


44.9 


29 f 


.34 


21.9 


6.7 


PS 


50 


365 


.976 


2.97 


62.6 


29 f 


.368 


28.8 


10 


a-PMMA 


43 


378 


.864 


2.35 


56.0 


29 f 


.33 


23.3 


7.7 


Pcarb 


51 


423 


.854 


2.81 


66.5 


73 


.22 


46.1 


5.7 


PaMS c 


56 


440 


.958 


3.98 


81.0 


56 


.32 


52.4 


11 



a a dT g JdP = T g V Aa/ACp. 

b Partially crystalline sample. 

c All quantities derived from the same 



ample. 



d 55 percent Styrene for T g , V g and Aa, 43 percent AC p. 
e Measurements by J. J. Weeks reported in Ref. 16. 
f Average over different sources. 



equation of Shih and Flory [36]. Although this equation is 
derived from data at temperatures well above T g , its nearly 
linear response apparently allows valid extrapolation to much 
lower temperatures. The value V g a t = 7.9 X 10~ 4 cm'7(g-rv) 
is slightly less than the average, 8.7 X 10 -4 , of the others for 
this polymer in table I of ref. [26] which are obtained from 
different sources of data at higher temperatures not encom- 
passing T g . Also the extrapolated value of V g = 0.904 cm 3 /g 
above essentially coincides with 0.905 in ref. [26]. The 
revised thermal expansion values give the positive value of k' 
shown in table 7. 

4.6. Comparison of Experimental Transition Rates 

In tables 6 and 7 the assumptions that dT g /dP could be 
replaced by {dT/dP)^ or 7VAa/ACp, respectively, are em- 
ployed. In table 8 values of these quantities are compared for 
each polymer. A similar table was prepared by O'Reilly [13] 
in 1962 for glass-forming liquids not restricted to polymers. 
Values of A/</Aa are also included here for comparison; 
however, agreement with dTJdP is not expected since the 
validity of the inequality 

dTjdP < Ak/Acx 

appears to be quite strong and general. In most instances, 
agreement between dTJdP, (dT/dP)^, and TV a/ AC P seems to 



Table 8. Comparison of experimental transition rates a 







PVT 




Dynamic 


Thermal 


Polymer 
































Ref. 


dT g /dP 


Ak/A(x 


Ref. 


(W/WL 


Ref. 


TVAa/ACp 


PIB 








61 


25 


59, 69 


24.1 


NR 


39, 60 




26 b 


60 


24 


39, 70 


24.8 


PP 


40 


20 








40, 71 


23.8 


PVAc c 


1 


21.2 


48.2 


64 


22 


1, 16 d 


22.1 


PVC 


47 


13.5 


33.5 


66 


18 


47, 29 


21.9 


PVC 


48 


14 


46 


67 


16.5 


48, 29 


28.6 


PS 


50 


31.6 


74.2 


68 


32 


50, 29 


28.8 


PS C 


9 


31 








9 


30.5 


PS 


55 


32 


73 






55 


34 


«-PMMA 


43 


23.6 


54.5 


67 


24.5 


43, 29 


23.3 


«-PMMA c 


9 


22 








9 


32 


Pcarb 








13 


44 


51, 73 


46.1 


PaMS 


56 


39 


52 






56 


52 



a All units in K/kbar. 

b Ak = 1.2 X 10 -5 bar -1 determined from dynamic compressibility. 

c All given determinations on same sample. 

d Measurements by J. J. Weeks reported in Ref. [16]. 



be within experimental error. Small differences may be antic- 
ipated because the conditions under which these quantities 
are evaluated may be vastly different. 

According to these results the Prigogine-Defay ratio 



r = Ak AC P /[7V(Ao0 2 ] 
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is essentially unity for natural rubber and PaMS. Unfortu- 
nately, we do not have a PVT value of dTJdP for the former 
to test the validity of the Ehrenfest equations [eqs (6) and (7)]. 
With PaMS it would appear that although neither of the 
Ehrenfest equations is obeyed, the Prigogine-Defay ratio is 
still unity, which is an atypical result. This implies that 
dT\ I dP for volume and entropy are equal, but dT g /dP is 
distinct. However, since data in tabular, or even graphical 
form, are not included in ref. [56], evaluation of these 
numbers cannot be scrutinized. Poor agreement for PVC in 
all cases is apparent; however, this may be a result of sample 
differences including the degree of crystallinity which is 
difficult to control in this polymer. Also poor agreement is 
noted for a-PMMA of ref. [9], where both quantities are 
obtained from the same sample. On the other hand, the data 
in the row above on the same polymer reveal good agreement 
including that with dTJdP (PVT) of ref. [9]. In all cases 
agreement is very good for polystyrene. 

These results indicate that dT g /dP = {dT/dP)^ is a valid 
relation and dTJdP = TV Aa/ &CP seems to hold most of the 
time. The validity of the first may be argued on a qualitative 
phenomenological basis (see sec. 2.) The second relation is 
evaluated at atmospheric pressure only. There is no apparent 
reason to assume that the approximation will be as good at 
elevated pressures. 

In section 4.3 we mentioned that PVT data on PnBMA 
suggest that the first Ehrenfest equation [eq (6b)] is a good 
approximation for this polymer. This result is tantamount to 
essentially no densification. (See table 3.) Unfortunately, we 
have no heat capacity values for this polymer, which are 
needed to check the second Ehrenfest equation [eq (7)]. In 
section 3 we noted the possibility of using the "rule of 
constant ACV' [20] to estimate the heat capacity difference at 
T g . (For a comparison of experimental and "bead" values of 
AC p on polymers, see ref. [28].) For PnBMA the molecular 
weight of the polymeric repeat unit is 142.2 g/mol. Assigning 
one bead to each of the two carbon backbones, and one to the 
oxygen atom, we obtain a total of three beads, which for 11.3 
J/(mol-K-bead) gives AC P = 0.24 J/(g-K). Taking this value 
along with those for Aa and T g from tables 3 and 4, respec- 
tively, and V g = 0.946 cm 3 /g from ref. [43], we obtain TVAa/ 
AC p = 20 K/kbar, which is in good agreement with 20.4 in 
table 3. Thus both of the Ehrenfest equations appear to be 
fairly good approximations for this polymer, along with a 
corresponding Prigogine-Defay ratio of nearly unity. (The 
value 1.2 is obtained for PnMBA. The average value obtained 
from table 8, exclusive of NR and PaMS, which were treated 
separately, is 2.1.) These results imply that both the density 
and entropy of PnBMA are essentially independent of forma- 
tion pressure, at least at low pressures. 

4.7. Pressure Dependence of k' and Limitations 
Imposed by Chemical Instability 

The previous discussions in this paper pertain to the initial 
values of k' or at least at very low formation pressures. Most 
of these are either tangent values at P' = (atmospheric 
pressure) or secant values obtained fromP' = 1 kbar or less. 
There are data in the literature, however, which include 
densifications obtained at different formation pressures. 

There are two important physical considerations in optim- 
izing the procedure to obtain "permanent," densified glasses. 
The first and more obvious, is to select and maintain the 
temperature of depressurization at temperatures sufficiently 



below T g . It is clear that the ambient conditions must be such 
as to maintain structural relaxation times which are large in 
comparison to the desired "lifetime" of the glass. Accord- 
ingly, high T g substances are preferable for room temperature 
stability. The second is to choose T , the temperature of 
isothermal pressurization, large enough that the equilibrium 
melt is always maintained during pressurization. Stated alter- 
natively, the inequality 



To > UP') 



(19) 



must be approximately satisfied, as illustrated in figure 1. 
This condition implies that the effective time of the pressuri- 
zation process must be large in comparison to the structural 
relaxation time at the final pressure P' . In cases where T < 
TgiP') there will be a much smaller contribution to the 
densification process when the condition T = T g (P) is 
approximated and exceeded during pressurization. This situ- 
ation is revealed by a leveling off in the volume as illustrated 
schematically in figure 4, where volume changes are plotted 
with respect to formation pressure at different pressurization 
temperatures. The densification is expected to be independ- 
ent of T at low pressures, when ineq (19) is satisfied, as is 
revealed by the coalescence of these curves with decreasing 
P'. Such a coalescence is not expected, however, when T < 
T g as is illustrated by the data of Shishkin [54] on polysty- 
rene. In figure 4 the dashed line represents the extension of 
the envelope established from arbitrarily large values of T . 




Figure 4. Illustration of the influence of the pressurization temperature 
T on the densification process. 

The dashed line is the envelope approached at large temperatures. 



One of the better experimental examples which illustrates 
the behavior shown in figure 4 is provided by the data of 
Shishkin on PMMA and PS. Formation pressures up to four 
kbar are applied; but not all of the pressurization tempera- 
tures are above T^P'). At the lower pressures, k' increases 
withP' as is indicated by the increasing slopes of Shishkin's 
volume-formation pressure curves, and as is shown in figure 
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4. This is the opposite of the trend for the isothermal com- 
pressibility K, which decreases with increasing pressure (see 
for example refs. [12] and [74]). The data of Shishkin, as well 
as those of Kirnmel and Uhlmann [75] on PMMA, show that 
some densification is possible with pressurization tempera- 
tures below T 0, but the effect is diminished as the difference 
between these two temperatures is increased. Other examples 
illustrating the dependence of densification on formation 
pressure are refs. [6], [9] and [55] on polystyrene, and [56] on 
PctMS. Yourtee and Cooper [6] observe a very slight decrease 
in the densification rate with formation pressure for poly- 
styrene over a 6 kbar range. For the same polymer, Weitz and 
Wunderlich [9] find a much larger dependence with the same 
trend, where the density gradually becomes nearly constant 
at 4 kbar. These trends are contrary to the marked increase in 
the densification rate with formation pressure observed by 
Shishkin on PS and PMMA and Ichihara et al. [56] on PexMS. 
k' does not necessarily have to tend to zero for the volume to 
be non-negative at large formation pressures. Using our 
definition of k', the densified volume tends to zero at arbi- 
trarily large formation pressures when k' is a positive con- 
stant. 

According to the experiments of Weitz and Wunderlich on 
polystyrene, there is a monotonic increase in density at a 
decreasing rate which the density seems to level off at 4 kbar. 
Thus, beyond this point the formation pressure would have no 
effect on the densification process. On the other hand, with 
most of the other investigations mentioned above, including 
those on polystyrene, it would appear that chemical stability 
is the limiting factor. Whether the reaction rate constant of a 
given rate process increases or decreases with pressure de- 
pends upon the sign of its corresponding activation volume 
[76]. In most cases it is expected that the total activation 
volume will be positive with a corresponding increase in the 
effective decomposition temperature with increasing pres- 
sure. This behavior may be complicated, however, by the 
different temperature and pressure dependencies of the var- 
ious decomposition modes, and, possibly by the initiation of 
new ones at elevated pressures. 

The important consideration here is whether, or not, the 
decomposition temperature and T g -pressure curves come 
sufficiently close at any point to limit the densification proc- 
ess. For example, with polytetrafluoroethylene the increase of 
decomposition temperature with pressure is only about 3.5 K/ 
kbar [77]. Although this rate is small, the decomposition 
temperature is sufficiently remote from the observed phase 
transitions, since its initial (atmospheric pressure) value is 
about 700°C. In addition, the melting and decomposition 
curves diverge with increasing pressure over the experimen- 
tal range of 28 kbar, investigated so far. 



In cases where the decomposition temperature 7\/ in- 
creases with pressure, T should also be allowed to increase 
with pressure to optimize the densification. With polytetra- 
fluoroethylene this process would appear to continue without 
bound because of the observed divergence mentioned above. 
In instances where T a and T g converge or intersect at a finite 
pressure, the densification would be essentially limited b) 
the effective intersection temperature as illustrated in figure 
5. Except for the polymer mentioned above, pressure de- 
pendent pyrolysis data are apparently non-existent in the 
literature 4 . 

The results in table 9 summarize an attempt to estimate the 
optimum densification on a few polymers beyond which ther- 
mal decomposition would occur. In the absence of reliable 
pyrolysis data at elevated pressures, we will estimate opti- 
mum densification by commencing isothermal pressurization 
at Tq = T & Since in most cases T a is expected to increase 
with pressure, this procedure should underestimate the maxi- 
mum densification as illustrated by the lower value of P' max 
obtained by the dashed line path in figure 5. 



4 Preliminary work of this kind on 18 polymers, including PMMA, PS, Pcarb, and PVC, has been 
reported by H. Eyring and coworkers in Ref. |7H|. In all of these the decomposition temperature 

increases nionotonically with pressure. We have been unable to obtain a final report on this activity. 



Pressurization 
in Melt 




Figure 5. Schematic illustration of temperature-pressure history used to 
optimize the densification process before the onset of pyrolysis. 







Table 9. Estimation of maximum densification f 


r om pyrolysis data 






Polymer 


Ref. 


A// kj/mol 


T K 


k(T) %/min 


7\,K 


T g K 


dTg/dP 

K/kbar 


P'max kbar 


k' Mbar 1 


^ j// max 

% 


PIB 


79 


218 


593 


0.268 


558 


202 


25 


14 


1.6 


2 


PP 


79 


255 


623 


.069 


606 


244 


20 


18 


7.0 


13 


PMA 


80 


155 


558 


.270 


515 


282 


18 


13 


5.8 


8 


PVAc 


81, 82 


112 


497 


5.58 


409 a 


304 


21.2 


5 


9.3 


5 


PVC 


83 


126 


508 


0.43 


458 


349 


13.5 


8 


6.4 


5 


PS 


80 


218 


608 


.163 


577 


363 


30.0 


7 


8.1 


6 


a-PMMA 


80 


230 


583 


.250 


552 


378 


23.6 


7 


7.3 


5 


PaMS 


80 


243 


546 


.276 


519 


448 


39 


2 


4.9 


1 



a Stated to be unstable at temperatures above 463K in Ref. [81]. During sample preparation [1] slight discoloration was observed after heating overnight in 
a vacuum at 403 K. 
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In table 9 T a is taken arbitrarily at the value for which the 
initial reaction rate constant k = 1%/hr, applicable to the 
total degredation process. Assuming Arrhenius behavior T a 
may be calculated from the relation 

1/T d = \/T - (R/AH)€n[(l/60)/k(T)] 

where R = 8.314 J/(mol-K) and AH is the activation energy. 
T is arbitrarly chosen from the closest data point to T d, which 
in all cases, but one (PVAc), involves extrapolation. These 
decomposition temperatures correspond to those given in 
table 7 of ref. [34], except the latter apparently apply to k = 1 
percent/min and, accordingly, are larger. The ceiling tem- 
peratures in the same table, which apply to the propagation 
mode at equilibrium, are apparently not relevant to the 
densification process. T g , dTJdP, and K f are selected from 
previous tables in this paper. P' max , the pressure corre- 
sponding to the onset of pyrolysis at T = T a and — (AV/ 
K) max , the corresponding maximum densification, are ob- 
tained from the simple relations 



P' 

-(AV/V) 



max 
max 



[T d - TMVidTjdP) 



' 9 
I 
max 



From these results it appears thatP' max varies inversely with 
T g ; however, no trend is apparent for (AV /V) max . 



4.8. Dependence of Refractive Index on Densification 

A reliable estimate of the change of refractive index on 
densification should be obtained by means of the Lorentz- 
Lorenz equation, 



l)/(n 2 +2) =Kp, 



(20) 



where n is the index of refraction and p the density. K 
depends upon the polarizability, which is expected to be 
essentially independent of formation pressure, or alterna- 
tively, the density at constant temperature and pressure. The 
relative change of index of refraction with formation pressure, 

8' = (1/rc) (dn/dP) Tt p, 

is obtained explicitly by differentiation of eq (20), viz, 

8' = (\/6n 2 )(n 2 - 1 )(n 2 + 2)*'. (21) 

Table 10 presents the results for polymers for which values 
for no (sodium D line) are available from ref. [35] with k' 
selected from our tables. The no values are converted to those 
at T g by means of the temperature coefficients given in ref. 
[84]. As seen from the table, these corrections are insignifi- 
cant. Since all of these values range between 1.48 and 1.58, 
a very slight (10 percent) error will be incorporated in 8' by 
taking the function f(n) = (l/6n 2 ){n 2 - l)(n 2 +2) as a 
constant, as revealed by the table. Accordingly, in view of 
the large experimental uncertainties in k' (35 percent for 
polystyrene), the additional uncertainties obtained on replac- 
ing eq (21) by the approximation 

8' = 0.4k' 



are slight. The values of 8' in the table however, are calcu- 
lated from eq (21). We do not have any direct experimental 
data giving the dependence of the index of refraction on 
formation pressure. 



Table 10. 


Estimation of change in refractive index from 


densification rate 


Polymer 


T g °C 


n D (Tf 
(°C) 


dnoldT* 
o c _, 


flD(Tg) 


f(n D ) 


k' 
Mbar" 1 


8' 
Mbar -1 


PIB 


-71 


1.51 c 


-0.0003 


1.54 


0.421 


1.6 


0.7 


NR 


-69 


1.52(25) 


-0.00037 


1.55 


.428 


1.9 


.8 


PP 


-29 


1.49 c 


( d ) 


1.51 


.400 


4.1 


1.6 


SBR 


-16 


1.53 c 


( d ) 


1.54 


.421 


5.8 


2.4 


PMA 


9 


1.47(20) 


( d ) 


1.47 


.373 


5.8 


2.2 


PnBMA 


20 


1.48(25) 


( d ) 


1.48 


.380 


0.7 


0.3 


PVAc 


31 


1.48(20) 


-o.oooi 


1.48 


.380 


9.1 


3.5 


PEMA 


64 


1.48(25) 


o 


1.48 


.380 


8.9 


3.5 


PVC 


76 


1.55 c 


( e ) 


1.54 


.421 


7.6 


3.2 


PS 


90 


1.59 c 


-0.00013 


1.58 


.449 


8.2 


3.7 


a -PMMA 


105 


1.49(20) 


-0.00012 


1.48 


.380 


7.4 


2.8 


PCHMA 


107 


1.51(20) 


-0.00013 


1.50 


.394 


12.6 


5.0 


Pcarb 


143 


1.58 c 


( e ) 


1.57 


.442 


5.1 


2.3 



a Taken from Ref. [35]. 
b Taken from Ref. [84]. 
c Temperature taken as 25°C. 



d dno/dT taken as 
e dno/dT taken as 



-0.0003. 
0.0001. 



These evaluations have potential application in optimiza- 
tion or adjustment of the refractive indices of plastic lenses 
by appropriately setting the molding pressure. The values in 
the last column in table 10 give the relative percent changes 
(An /n) resulting from a moding pressure of 10 kbar. For PS 
and PMMA, which are common constituents for plastic len- 
ses, n would change by 4 and 3 percent, respectively. 
However, it was estimated in the last section that thermal 
decomposition of these polymers would limit the pressuriza- 
tion to 7 and 8 kbar, respectively. In these analyses isother- 
mal pressurization is considered at the decomposition tem- 
perature. If this temperature increases with pressure as indi- 
cated by ref. [78] for PS and PMMA, an additional increase 
in their refractive indices could be obtained by appropriately 
increasing the temperature during pressurization. 

5. Conclusion 

Several methods have been evaluated to estimate the den- 
sification rate, k\ applicable to glass formation by isobaric 
cooling at constant rate. Other than the direct measurement 
of the volume difference in the glass, k' is always computed 
from an expression involving the difference between two 
transition rates, dT g /dP and dT^/dP. The hole theory is 
shown to be sufficiently accurate in estimating dT^/dP for 
the 23 polymers evaluated except for possibly those of di- 
methyl siloxane and a-methyl styrene. With these it is not 
clear whether the discrepancies result from experimental 
error or lack of generality in the application of the theory. 
Although dTg/dP is only evaluated experimentally, there are 
independent alternatives. The simplest of these involves the 
differences between thermal expansions and heat capacities 
at T g for liquid and glass at atmospheric pressure only. 

The principal problem in the estimation of densification 
using these procedures appears to be the large amount of 
experimental uncertainty in all of the relevant quantities, in 
particular, the compressibility. Since the expression for K r 
involves the difference between two quantities of similar 
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magnitude, even small experimental errors may have a pro- 
nounced effect. Accordingly, it is difficult to assess the 
relative merits of the different methods employed here, in- 
cluding the application of the hole theory. 

The results of these analyses appear to have 1 practical 
applications. Densifying glasses produces a hardening effect 
as revealed hy an increase in moduli. However, these effects 
do not appear to be as pronounced, in particular viscosity or 
relaxation time, as those obtained at the same volume, tem- 
perature, and pressure in the glass by commensurately de- 
creasing the cooling rate at atmospheric pressure. This pro- 
cedure however is usually not practical because of the large 
times required for glass formation. According to one investi- 
gation it is possible to optimize the ultimate properties 
through the appropriate adjustment of the formation or mold- 
ing pressure. More work is necessary to establish the general- 
ity of this result and to determine the formation pressures for 
maximum yield stress. Moreover, the relation between the 
refractive index and densification quantity presented could 
be used to quantitatively regulate the refractive index of 
lenses through appropriate adjustment of the molding pres- 
sure. The maximum value would appear to be limited by 
chemical instability at the high temperatures necessary to 
exceed T g (P), which increases with pressure. The simple 
relation given does not include the influence of densification 
on optical dispersion. Again, experimental work is required 
to assess the validity of our estimates and the possible 
influence of densification on dispersion. 



This work was supported in part by the National Science 
Foundation under Grant DMR 75-15401. 



values were chosen so as to magnify the effect. For the glass, 
the approximation will be even more adequate. With a pres- 
sure difference of 2 kbar and (k) = 4 X ]0~ 5 bar -1 , the 
magnitude of the relative volume change is the same and 
identical conclusions are obtained as for the a-term. 

Denoting the initial temperature in the melt as T and the 
final temperature in the glass as T '/, with the initial pressure 
taken as zero, we obtain instead of eq (3), when the pressure 
dependence of the a's and the temperature dependence of the 
k's are taken into account, the following expression: 

(Va - V B )/V = f '[><i(To, P) ~ x g< AT,, P)] dP 

- f T %l«i(T, 0) - a„, c (T, P')] dT + f T \p.l<xi{T, P') 

1 g a 

- aAT. 0)] dT - j'; f ' m [a g , b (T, 0) - a 9 ,-(7', P')] dT 

(A-l) 

where the subscripts b and c pertain to the low and high 
pressure glasses respectively. This choice conforms with the 
nomenclature used in refs. [1] and [14]. In the linearized 
derivation, a and K for both liquid and glass are taken to be 
constants. Accordingly, there is no distinction between the 
values of (X g £ and <x g , c , and, similarly, Kg^ and Kg tC . 

To proceed further, we make use of an explicit equation of 
state. It is most convenient to employ the extensively tested 
Tait relations for both melt and glass. To recapitulate the 
pertinent equation (14), [50]: 



6. Appendix 

As indicated earlier, eqs (5) in the text involve a lineariza- 
tion of pertinent quantities and proximity to the transition line 
T g (P) in figure 3. The experimental data found in the litera- 
ture often do not satisfy these conditions. Hence we reconsi- 
der here the processes depicted in figure 3, by replacing the 
simplifications adopted earlier by the more general form. 
This will not only permit the prediction of densification 
effects under more extreme conditions, but also allow us to 
gage the quantitative validity of the linearization. Clearly, a 
knowledge of the equations of state is required for an explicit 
evaluation, but is certainly not available for the wide range of 
systems discussed here. However we shall be able to present 
typical numerical illustrations using PVAc where appropriate 
data exist [1, 14]. 

From the definitions of the coefficients a and K it follows 
that: 



k(P, T) = C{[P + B] X [1 - C€n(l + P/B]}- 
e 

B(T) = aexp(-bT) 



a(T, P) - a(T, 0) = Pk{T, P) d€nB/dT = -bPK. 



(A-2) 



The last transforms eq (A-l) into 



(Va ~ V B )/V C = fo'MTo, P) ~ K„, c (T f , P)] dP 

- ffgfoatf, 0) - a g , c (T, 0)] dT - /Ma^r, 0) 



V/V = exp U T To <*dT]; V/V = exp [-£ KdP] 

where the subscript indicates an initial value. However, it 
will be accurate enough to omit second and higher powers in 
the expansions of the exponentials. Considering an average, 
temperature independent value for the liquid, (a) = 8 X 
10~ 4 K _1 , and a temperature interval of 100 degrees, we 
obtain for the integral a value of 8 X 10" 2 . Thus the 
quadratic term changes the total result by 0.34 percent. The 



cCr.ofldr-wjI* «&,?) 



s(P') 



XdT 



bo,cP' r« p '>K atC {T,P')dT, (A-3) 



where the first two terms will predominate. 

We now proceed to evaluate the integrals in eq (A-3) which 
are identified as follows: 



295 



h = /o fain, p) - K g ,c(T f , p)]dP^\KP' 



References 



h = SffiMT, 0) - a 9 , c (T, 0)]dT 

1 9 



ku(dT g /dP)P' 






/;; [a g , b (T,0)-a g , c (T,0)]dT- 



-btP'f'^^AT^P^dT^O 

l 9 

-b B J"Sl^K a , c (T,P')dT^0 



(A-4) 



where the terms to the right of the arrows are the correspond- 
ing linear approximations used in eqs (5) in the text. 

Since C = 0.0894, the compressibility may be written in 
good approximation, and consistent with the expansion of the 
exponentials above, as C/(B + P). Thus we find for/j, 



I x = C£n 



1 + P'/BiJTq) 
\+P'/B , c (T f ))' 



1 2 and I 3 are evaluated by expressing the atmospheric pres- 
sure volume as V = A + B^T + C qT 2 for which a — (B + 
2CoT)/V av is a good approximation. V av is taken as the 
average of the two bounds. Finally the general integral corre- 
sponding to 1 4 and I 5 is 

J k(T, P) dT = (C/P) {T 

+ (1/6) €n [P + a exp(-bT)]} + f(P). 

Using the parameters for PVAc given in tables 1 and 2 
of ref. [14], the values of the integrals and their linearized 
counterparts are summarized in table A-l. From ref. [1] Tg(0) 
= 30.7 °C and T g (P') = 48.0 °C. T and T/are taken to be 
90 and °C. These two temperatures are considered to be 
sufficiently remote from T g (P), to be characteristic of the 
equilibrium and glassy states, respectively. The total relative 
volume differences are given at the bottom of the table 
followed by the corresponding densification rates. The differ- 
ence between the two values of k' amounts to about 4 
percent, which is quite satisfactory, since the experimental 
error on this quantity appears to be considerably larger. This 
good agreement seen in table A-l arises however from a 
cancellation of approximation errors. Moreover, it is gratify- 
ing, that the value k' = 8.8 Mbar -1 , based on the Tait 
equation is essentially identical to that obtained by directly 
measuring the volume difference. (See sec. 4.3.) This illus- 
trates the satisfactory performance of the analytical expres- 
sions in representing the experimental data. 



TABLE A-l. Values of integrals in eqs (A-4) 



i 


Integral Value, /j 


Linear Counterpart 


1 
2 
3 
4 
5 


2.042 X 10-2 
-0.649 X 10-2 

0.106 X 10-2 
-0.552 X 10~ 2 
-0.242 X 10-2 


1.465 X to -2 
-0.733 X 10"2 





(Va-VbWo 
k' (Mbar *) 


0.705 X 10- 2 
8.8 


0.732 X 10- 2 
9.2 
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